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Abstract 

A description of the DEPOSIT computer code is presented. The code is intended to calculate total and m-fold electron-loss cross 
sections (m is the number of ionized electrons) and the energy T{b) deposited to the projectile (positive or negative ion) during a 
collision with a neutral atom at low and intermediate collision energies as a function of the impact parameter b. The deposited 
energy is calculated as a 3D-integral over the projectile coordinate space in the classical energy-deposition model. Examples of the 
calculated deposited energies, ionization probabilities and electron-loss cross sections are given as well as the description of the 
input and output data. 
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1. Program summary 



^ Program title: DEPOSIT 

Licensmg provisions: GNU General Public License version 3 
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■No. of lines in distributed program, including test data, etc. 



Distribution format: tai-.gz 

^ Programming language: C++ 

Computer: Any computer that can run C++ compiler 

Operating system: Any operating system that can run C++ 
. compiler. 

o ■ 

(N Classification: 2.4,2.6,4.10,4.11 
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Nature of problem: For a given impact parameter b to calcu- 
late the deposited energy T(b) as a 3D integral over a co- 
ordinate space, and ionization probabilities Pm{b)- 

For a given energy to calculate the total and w-fold 
electron-loss cross sections using T{b) values. 

Solution method: Direct calculation of the 3D-integral T(b). 
One-dimension quadrature formula of the highest accu- 
racy based upon the nodes of the Yacobi polynomials for 
the cos^ - X e [-1,1] angular variable is applied. The 
Simpson rule for the tp e [0, In] angular variable is used. 
Newton-Cotes pattern of the seventh order embedded into 



every segment of the logarithmic grid for the radial vari- 
able r G [0, oo] is applied. Clamped cubic spline interpo- 
lation is done for the integrand of the T(b). 

Bisection method and further parabolic interpolation is ap- 
plied for the solving of the nonlinear equation for the to- 
tal cross-section. The Simpson rule for the m-fold cross- 
section calculation is applied. 

Running time: For a given energy, the total and m-fold cross 
sections are calculated within about 15 minutes on 8-core 
system. The running time is directly proportional to the 
number of cores. 

2. Introduction. 
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Electron loss (also called projectile ionization or stripping) 
of positive and negative ions, colliding with neutral atoms, is an 
important collision process playing a key role in many problems 
of plasma physics |[ll]-|S] and physics of accelerators itrll- lfioll . 
For example, electron-loss (EL) cross sections are required for 
estimation of the ion-beam lifetimes and vacuum conditions 
(residual-gas pressure and concentrations) in accelerators and 
storage rings, for energy losses of high-temperature plasmas 
and for heavy-ion beam-probe (HIBP) diagnostics of plasmas in 
tokamaks and stellarators. General properties of various atomic 
processes occurring in the collisions of heavy ions with neu- 
tral atomic and molecular targets are investigated in detail and 
presented in various review articles and books lfll1l - ll20h . 

At low and intermediate energies, a multiple-electron losses 
contribution to the EL cross sections can be very large [21]- [231] 
(up to 50% and more), but with the energy increasing a contri- 
bution of multiple EL processes goes down and single-electron 
losses become dominate. At high and relativistic energies the 
EL cross sections can be calculated in the relativistic Born ap- 
proximation II24II - II27II . At low and intermediate energies, the 
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total EL cross sections can be calculated in the classical ap- 
proximation using the CTMC approach ll22ll . ll28ll - ll29ll (Clas- 
sical Trajectory Monte Carlo) or the energy-deposition model 

JH-il. 

In the CTMC approach all projectile and target electrons are 
considered as classical particles described by a system of a few 
hundred coupled Newton-type first-order differential equations 
with a few thousand trajectories. The corresponding computer 
code is quite complicated in use and takes a very long run-time. 

The DEPOSIT program was created as an alternative code 
to the CTMC one and is intended to calculate the energy T{b) 
deposited as a function of impact parameter b to the projectile 
(positive or negative ion) during a collision with a neutral atom. 
The obtained T(b) values are then used for total and m-fold EL 
cross-section calculations (see Section [3]l. The code is based 
on the classical energy-deposition model suggested by N. Bohr 
1 34] and developed later in ll35ll - ll36ll for multiple-electron ion- 
ization of atoms by positive ions. In the DEPOSIT program, the 
problem of obtaining the energy deposited to the ions by neutral 
atoms is reduced to the calculation of a 3D-integral over the ion 
coordinate space. The code requires much less computer time 
and provides approximately the same accuracy of the EL cross 
sections as the CTMC method, i.e. a factor of 2. 

In this paper we describe a numerical implementation of the 
energy-deposition model in the DEPOSIT computer code and 
give a brief user guide. One important property of the code is 
that it allows one to simulate electron-loss cross sections for 
an arbitrary projectile (neutral atom, positive or negative ion) 
colliding with arbitrary atomic target. The validity conditions 
of the code are also discussed. 

The article is organized as follows. In Section |3] we briefly 
describe basic ideas and physical points of the developed 
method. In Section |4] we use these results and explain the nu- 
merical implementation of the model. Some basic examples on 
how to use the DEPOSIT code from the user's point of view are 
set out in Section|5] Finally, conclusions are given in Section|6l 
Atomic units are used throughout the paper. 



3. Physical model 

The energy T(b) transferred to the projectile ion by a neutral 
target atom as a function of the impact parameter b is given by 
a 3D-integral over the projectile coordinate space 113 311 : 



T(b)^Yj J Pyi'-)^Ey{p)dh. (1) 



In the sum index y labels principal n and orbital I quantum num- 
bers of atomic shell, y - nl. The summation is over all occupied 
shells of the projectile ion. 

In the DEPOSIT code, the projectile electron density Py{r) 
for the y-shell is presented by the nodeless Slater wave function 
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r{2i^ + 1) 



1/2 



r^e-f^', (2) 



and normalized to the total number of the projectile electrons 

(3) 



Here T(x) is the Gamma function, /? and /i are Slater parame- 
ters il, ii. 

In the classical approximation the kinetic energy, deposited 
to the shell y by the target atom can be expressed in the form: 



{Uy + t^Myf- Uy 



(4) 



where Uy is the orbital velocity of the electron in the y-shell and 
Auy denotes the energy gain by the electron. The values Uy and 
Auy are obtained from the classical equations 



J-co dp 



dU(R) dR 



dt. 



(5) 



(6) 



dp J_oo dR dp 

where ly is the binding energy of the projectile y-shell. The 
function U (R) denotes the field of a neutral atom at a distance R 
from its nucleus and is calculated as a sum of three Yukawa- 
type potentials with five fitting parameters A,, a, obtained from 
the Dirac-Hartree-Fock-Slater calculations [,39,1 : 
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2 -1- uj. 



(7) 



(8) 



Here v denotes the velocity of the projectile (the straight-line 
trajectory approximation is used) and is the relative velocity. 

The impact parameter p is a function of b and r and in spher- 
ical coordinates (r, 0, (p) has the following dependence in the 
frame of the moving projectile: 



- (b- r cos 6) + (r cos (f sin 6) 



(9) 



The energy gain Auy in eq. ^ diverges at small impact pa- 
rameters p and must be renormalized 113 311 . After normalization 
of Auy the deposited energy AEy is finally determined depend- 
ing on the relation between ion velocity v and orbital velocity 



AEy(p) = AE;(p)e(Uy -V) + AE;(p)e(V - Uy), 



where 



AE<(p) : 



AE^ip) 
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The step theta-function 6(x) is defined in a standard way 
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In order to avoid the thresholds in every point v - Uj, due to the 
discrete behaviour of the 6(x) (i.e. to make the deposited energy 
AEyip) a smooth function), the step theta-function is replaced 
by the right-hand side Fermi-smearing function 



n/(x) 



1 



c-kx 
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(14) 



The n/(x) becomes the 9{x) when ^ — > oo. Parameter value 
^ = 3 is taken by default. It is an input parameter and can be 
changed, if it is needed. Thus, 

hEyip) = ^E;{p)nf{uy -v) + AE^{p)nf(v-Uy), (15) 

(v) 

The constant variable A^^^ depends on the mean radius ry of 
the projectile shell and on the atomic radius Ra of the target 
atom 



(y) 
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ry < Ra, 
ry > Ra, 



(16) 
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where ry - pyly ' . The mean radii Ra for neutral atoms are 
given in ll40ll . The validity conditions are discussed in |31]. 

The function F{y) in eqs. ( fTTT i and ( fTSl i is a result of using 
the potential U (R) in the form O and is given by 



Jo 

It can also be represented as 



(y2 + x2)3/2 y2 + 



F(y)=yKiiy). 



dx. (17) 



(18) 



where Ki{y) is the modified Bessel function of the second 
kind (see [41], p. 375, eq. 9.6.2, [42]). Integral ([T8]l does con- 
verge for all positive y. It has the following expansion for y 
close to zero (which can be obtained from ref. [Ull, eq. 9.8.7, 
p. 379 and eq. 9.8.3, p. 378) 

F{y ^ 0+) = 1 -H In I -H 0.03860786j/ + 0(/). (19) 

To speed up the calculations it is interpolated by using of 
clumped cubic splines ll43ll . p. 97. 

Once the T{b) energy is determined, the total EL cross- 
section can be obtained directly as 



O-tot 



(v) = 2< 



nbi 



m— 1 



where ^^max is (numerically) found from the equation 
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and /] denotes the (first) ionization potential of the projectile. 
The m-fold EL cross-section cr,„(v) can be calculated as 
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o-m(v)-2n P,n(b)bdb, 



(22) 



where Pm(b) is a probability to ionize m electrons from the pro- 
jectile. 

In the CTMC method the probabilities Pm{b) are calculated 
numerically. To calculate P„,{b) in the present approach the 
Russek-Meli model |44] is used. In this model the Pmib) prob- 
abilities are given in a close analytical form as functions of the 
deposited energy T{b) and the ionization potentials 4 of the 
projectile 

Pmib) = i^^S,Mb))l g i^. j 5,(77,(fo)), (23) 
S„{x) = 2l<"'-'>^2l7r''"/21x<3'"-2)/2/(3,^^ _ 2)11 (24) 



(25) 



Here denotes the binomial coefficient, {a} is the integer part 
of a. Double factorial means m\\ - m(m - 2){m - 4) . . . p with 
p = I for odd m and p - 2 for even m correspondingly, and 
Sm(x < 0) = by definition. 

The probabilities Pmib) satisfy the normaUzation condition 
at every impact parameter b 



YjP'nib) 



(26) 



m=l 



as it follows from the equations (|23]|-(|25]|. The variable A^' here 
means the number of active electrons considered to be ionized. 
Actually, A^' in equations ( l23T l, ( l26b should be equal to the num- 
ber of projectile electrons A^, but due to the rapid convergence 
of the sum ( |26] | it can be cut off at A^' < A^ value. 

Using the deposited energies Tib) calculated in the DE- 
POSIT code and the ionization potentials h of the projec- 
tile ion, the m-fold EL cross sections (T,„(v) can be calculated 
straightforwardly by applying equation (l22t . 

4. Numerical procedure 

From the numerical point of view the challenge for the evalu- 
ation of the deposited energy ([TJ is a three-dimensional integral 
over the projectile-ion coordinate space. For a given y-shell it 
can be written in general form as 

Tyib) =1 dr I dx I d<pRir)Fib, r, x, (p). (27) 
Jo J-i Jo 

For simplicity, y-notation is skipped in the right-hand side of 
the integral. Parameter b is fixed, and functions R{r) and 
Fib, r, X, (fi) are smooth, real and do not have any singulari- 
ties. Function Fib, r, x, tp) is nothing but AEyip) given by the 
expressions (fTTT l and (fT2l l. and the radial function R{r) has a 
Slater-type behaviour and represents the electron-density con- 
tribution dJJi 

Rir)^NiC\r^^'e-^'^'. (28) 

Here C\ is a normalization factor from the equation ([2), A^i is 
the number of electrons in given y-shell. 



In the expression (IZTT i all three integrals over the variables r, x 
and if have a constant upper and lower limits. It means that each 
can be split out and considered separately like one-dimensional 
integral with a given set of fixed parameters, namely 



,rx(b, r,x)- I d(pF(b, r, x, ip), 
Jo 

- J dxFhr.xib, r, x), 

f 
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Ftr{b,r) 

Ty(b)= I drR{r)Fbr{b,r) 



(29) 

(30) 
(31) 



4.1. Integration over the angles 

In the integral ( |29] | the dependency on (f is expressed only 
by cos^^ as it follows from the equation (|9]l. Due to the sym- 
met ry of the cos^i^ function, implementation of the Simpson's 
rule 14511 . ll46ll to the integral (|29t can be done in the following 
way 



/(cosV)# ~ ^ + /o + 4/i + fy) , 



(32) 



NJ2 



S^ = Yj (4/2^+1 + 2/2*) > fk = /(cos2 ^), (33) 



2Nc 



where Nc must be even and /(cos^ <f) - F(b, r, x, (p). Discrep- 
ancy estimation for eq. (l32T i is given in ref. ll46ll . 
Integral (l30l l can be done by the following rule 
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(n + l)(n -1- 2) 



jAt/(^i), (34) 



k=i 



8(n+ 1) 



(« + 2)(l-x2)2[(Pi''"(x,))']2 



(35) 



with an error estimation given in ref. [46]. Here, P^^'^\x) is 
a Yacobi polynomial, x^ is the A:-th root of the polynomial 
P^^'^\x), fix) = Fbrxib, r, x). For the choice of the « = 53 nu- 
meric error is equal to A- 10 where numeric factor A depends 
on a given system, i.e. the function Fhrxib, r, x). 

To speed up the computing procedure based on formulas 
( l34b . ( l35T l preliminary calculations of the Yacobi polynomial 
roots Xk and coefficients have been done. These two data 
sets have been directly built in the present code. 

4.2. Integration over the radius 

Due to the Slater-type behaviour of the radial function (l2Ft . 
integrand R{r)Fi,r{b, r) decreases exponentially. The main con- 
tribution to the integral (|311 function R{r)Fi,r{b, r) makes in the 
range [0, r^ax] ■ The infinity limit can be cut off in the rmax point. 
The value of r,nax is an input parameter, which is ordinary taken 
approximately equal to 100, depending on the projectile ion 
size. Input parameters are described in Section |5] 

To do integral (ISTT i the logarithmic grid with an input param- 
eter Ngrid is used 



r. - r Ji-Nsrld)lto 
' I — ' max^ 5 



< / < M 



grid-> 



to = 32, (36) 



and a Newton-Cotes quadrature formula for the seven-points 
interval 

'''' = ^ ''^^^"'--P ^ 27/(x,,.)+ ^^^^ 

272/(x,^3 ) + 27/(x,^4) + 216/(x,.^5) + 41/(x,-+i)], 



is applied for every range [r,, r,+i] (see ref. 114 111 , p. 886, 
eq. 25.4.16). Here x, - r, and points in the interval [x,, x,+i] 
are distributed uniformly. 

In this way integral (ISTT l can be done by the following rule 



Jo 



'f(r)dr^ 2 /r'. 

1=0 



(38) 



4.3. The total cross-section calculations 

The total-cross section cr,o, is calculated from equation (|20] |. 
where the /^max parameter is taken as a numerical solution of 
equation ( |2TI ). Function T{b) monotonically decreases with in- 
creasing of parameter b near the /j^ax point (a typical example 
of T{b) behaviour is given in Figure [TJ- 

To solve equation (1211 1 under these conditions the bisection 
method is applied within a given range [bi,b2], where bi and 
b2 are input parameters. The bisection procedure is stopped 
when the condition 



^2 



• bf < s. 



10 



-3 



(39) 



is satisfied, where / is the number of iteration. Then the 
quadratic interpolation procedure is used to find the local mini- 
mum of the (r(/7)-/i)2 function near the b^aRx value using three 
abscissa points 



■s/4. 



bo + s/4, 



(40) 



where bo - (b^^"^ + b2°^)/2, io is the final bisection step con- 
trolled by the condition (39[ . Once the b^^^ parameter is de- 
termined, the cr,o, value is calculated by applying the equa- 
tion ( l20l i. 



4.4. The m-fold cross-section calculations 

The m-fold cross sections cr,„ are calculated by formula ( 
where the integration is done from up to femax- To do inte- 
gral ( l22l i one, first of all, needs to calculate T(b) values in every 
point of the integration grid as it follows from equations (l23l l 
and (l25T l. 

The grid is taken with a fixed step h < 0.01 on the [0, femax] 
range. The number of points 

Np = [(1 + {^max/O.Ol)) div 2] ■ 2, (41) 

and the grid points are 

xj^jh, j = 0,1... Np, h^b^,JNp. (42) 

The integer part of a in equation (1411 1 is labeled as {a}, and the 
integer division operation is noted as div. 
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The Simpson rule 



f(x)dx « - [/(xo) + F2 + F4+ f(xN^)] , (43) 



^2 = 2 ^ f(X2j), F4=4^/(X2y-,) (44) 

is applied to do integration (l22l i on grid (l42t . 

The probabilities P„,{b) are implemented in correspondence 
to equations (l23T l - (l25T l with the following boundary conditions 



^'m(^ < eo) 



P„l{b > Mao) = 



1, m — I, 

0, 1 < m < A?', 



0, 1 <m< N', 



(45) 



(46) 



1, m = N', 

10^*', A^' is a cutoff value for the 



where eo = 10" and Mo 
sum ( |26] |. 

To avoid divergence in the numeric calculations of factorials 
and powers of 2, n and x in equations ( l23T l and (l24l i the follow- 
ing approach to Pm(^) calculation have been implemented 



Pm(b) 



N' 



n-1 



1+ ^ exp(Af(fe)) 



/= 1 Ji^m 



(47) 



A™(^) = S f + S f + S f + S f + S fib) (48) 



N'-m 



/;i=A^'-/+l ^:2=m+l 
m A'''—/ 



A:i=/+1 k2=N'-m+l 
3m-2 3!-2 

5f = ^ In^i - ^ In ^2, 

k,+=2 fc+=2 



(49) 



(50) 



(51) 



Sf = ({(/- l)/2) - {(m - l)/2))ln2, (52) 



5r = ({//2)-{m/2))ln7r, 



(53) 



Sfib) = (3il2 - l)\nrji{b) - Omjl - \)\xyn^{b). (54) 

Summation in the first sum of equation (ISTl i over k\ is done 
for all even positive k\ < 3m - 2, if m is even, and for all odd 
positive fei, if m is odd, as well as summation over even/odd 
positive k2 for even/odd / in the second sum of equation ( BTT l 
correspondingly. 

Integral (l22l i for every cr,„ is done on the same grid (l42l l 
[0, ^max]- An example of the P,„{b) probabiUties is given in Fig- 
ure |2] 



4.5. Parallelization 

Integral ([TJ representing the deposited energy T(b) after 
all implementation steps is finally reduced to the summation 
rule (l38T l. It can be broken down into the Nmre subsumes where 
Ncore is the number of MPI (Message Passing Interface) threads 
(cores). For parallel version of the present code the OpenMPI 
library is used [i47ll . 

In the parallel version integration in (l3Ft for a given ^-th 
core, < ^ < Ncore, is done by the rule 



hlil-i 

'■='1 [?] 



where 



h[q\ 



^ core ) 



\ii{q+\], ii{q+\]<Ngrid, 



\Norid, Otherwise. 



(55) 



(56) 



(57) 



The integer part of a is labeled as [a]. After the partial summa- 
tions ( l55T l on every core the final result 



q=Q 1=0 



(58) 



is collected by calling the Allreduce MPI function. 

As one can see the parallelization having done in the 
way ( fSST i does not harm the rule (l37t . This remarkable prop- 
erty combines high accuracy with high performance. 

5. Input/output 

5.7. Input 

All input parameters for the program are set up by means of 
one input text file with the proper keywords. Every keyword 
can be put in arbitrary place of the input file without a definite 
compliance of the order. The program is case-sensitive to the 
keywords. Comments can be added in CjC++ programming 
language style. 

The keywords are: Vi - velocity of the projectile in frame 
of the neutral atomic target given in atomic units, ZA - atomic 
charge and RA - effective radius of the atomic target, see equa- 
tions ( fTTT ). ( fT2l i and ( fT6b . Each of these parameters should be 
followed by one or several spaces and a real number (the value 
of the parameter). For example, for the BcP'^ + O collision at 
the energy E - 2.5 MeV/u they are 



Vi 
Za 
Ra 



10.0 
8.0 

1.239652 



Two values A\ and A2 for the equation O are given after the 
keyword A_exp. Three values a\, a2, aj, for the equation Q 
are given after the keyword alf _exp. 



A_exp 
alf _exp 



0.0625 
14.823 



0.9375 
2 . 0403 



0.0 
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First integer value after the keyword shells sets up the num- 
ber of the ion y-shells. The number of electrons Ny, the values 
of the coefficients C\y, ^y, j3y for the density (l28l l and binding 
energy ly must be put in as well. For every y-shell it is more 
convenient to set up the five parameters from a new line. Rules 
for calculation of the Ciy, iJ.y,/3y parameters are given in the ref- 
erences liSlll . 1 381. The binding energies ly are p ut in electron- 
volts (eV). The values of ly can be found in H, ll48ll-||52|. 

Shells 7 



//N 


CI 


mu 


beta 


KeV) 




8 


7 


778 


4 





2 


5625 


33 


11 


// 5sp8 


10 


121 


07 


3 


5 


4 


8143 


106 


044 


// 4dl0 


8 


956 


5335 


3 


5 


8 


0714 


215 


52 


// 4sp8 


10 


2252 


830 


3 





11 


6167 


801 


35 


// 3dl0 


8 


5404 


971 


3 





14 


9167 


1152 


6 


// 3sp8 


8 


3951 


539 


2 





25 


925 


5542 


13 


// 2sp8 


2 


831 


405 


1 





55 


7 


37455 


41 


// ls2 



The keyword b_range gives bmin, ^max, and Ab for calculat- 
ing and plotting T{b). 

II T(b) curve 

b_range 0.0 3.0 0.01 

The step parameter Kb does not need to be the range devided 
by an integer; /7-points are chosen by the rule bi - b„ii„ + iAb, 
i = 0,1,2. ..k,bk < b,„ax, bk+i > b,„ax. 

To calculate the total cross section cr,^, one needs to give 
a keyword Sigma_tot with three parameter values: bisection 
start range [bi,b2] (see Section l43T l and the first ionization po- 
tential /i [eV]. 



// Total cross section 
Sigma_tot 0.0 11.0 



34.45 



To calculate m-fold cross sections cr,„, where I < m < N', 
one needs to give a keyword Sigma_m_f old with A^' + 1 val- 
ues: the value of A^' parameter and A^' values of m-th ionization 
potentials /„, [eV] (see equations (I22l i-(|26]| and Section l4!4] i. 

// m-fold cross section 
Sigma_m_f old 30 

34.45 48.40 62.35 76.30 92.53 

107.1 139.2 155.7 232.2 266.0 

299.8 333.6 367.4 401.2 437.8 
472.0 506.2 540.5 672.6 710.3 

747.9 785.6 839.2 878.2 975.9 
1016.0 1586.0 1688.0 1790.0 1892.0 

One important point is that in order to calculate m-fold cross 
sections one must also calculate the b^-dx value from equa- 
tion {2TI . It is needed to cut off the integration grid in equa- 
tion (l22T l. Thus, the Sigma_m_f old keyword must be given in 
the text file with the keyword Sigma_tot and all their data. 

The remaining parameters should be kept with the default 
values. 



rgrid 

ksmear 

cosN 



70.0 600 
3.0 

54 



30 



Parameter rgrid sets up the r^ax, Ngrid and to values for equa- 
tion ( l36l l. The r,nax corresponds to the physical infinity for the 
atomic radial variable. Parameter ksmear fixes the k value in 
equation (fT4l l and parameter cosN sets up the number of points 
Nc in equations (l32ll- (l33l l. Value of the cosN parameter must be 
only even. 

5.2. Output 

The name of the input file is passed to the program input with 
the name of the output file via command line when the program 
starts. The output of the program is sent to console. 

> deposit f ile_namel_inp f ile_name2_out 

The program writes a text file with the deposited energy 
called energy_Tb_[f ile_najiie2_out] , if the b_range key- 
word is given. The file consists of several columns. The first 
column contains the /7-values, the second column is the T{b) 
energy and the rest of the columns are the deposited energies 
for every y-shell. 

The program also writes a text file with the probabiU- 
ties Pmib) called probability_Pm_[f ile_name2_out] , if the 
Sigma_m_f old keyword is given. The file consists of several 
columns. The first column is the /7-values in eq. (l23T l and the 
other ones are values of the P„,ib) probabilities for the corre- 
sponding b points and 1 < m < N'. 

All input parameters are printed out to the console (for the 
user's control) and the program starts the calculations. 

**** DEPOSIT CODE (VERSION 1.27/P) **** 
Number of MPI Threads = 8 
Input reading from file: ba2+_o.txt 
Vi = 10.0 Za = 8 Ra = 1.23965 



A_exp : 
alf _exp : 



0.06250 
14.82300 



0.93750 
2.04030 



Shells 
N 

8.0 
10.0 

8.0 
10.0 

8.0 

8.0 

2.0 



= 7 

CI 
7.778 
121.070 
956.534 
2252.830 
5404.971 
3951.539 
831.405 



mu 
4.0 
3.5 
3.5 
3.0 
3.0 
2.0 
1.0 



beta 
2.562 
4.814 
8.071 
11.617 
14.917 
25.925 
55.700 



0.00000 
0.00000 



I[eV] 

33.11 
106.04 
215.52 
801.35 
1152.60 
5542.13 
37455.41 



I[au] 




u 




Neff 


1.22 


1 


56 






3.90 


2 


79 






7.92 


3 


98 






29.45 


7 


67 






42.36 


9 


20 






203.67 


20 


18 





0511 


1376.46 


52 


47 





0019 
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2.317627 0.002461 

(r_max = 70.0 Ngr = 600 tO = 30.0 2.318970 -0.001648 

ksmear =3.0 cosN = 54) 2.318298 0.000405 



Slater w.f . normalization test: 



Interpolate : 



P. 


.1 


1 


0000 


2.318384 


0.000143 


p. 


.2 


1 


0000 


2.318634 


-0.000622 


p. 


.3 


1 


0000 


2.318884 


-0.001386 


p. 


.4 


1 


0000 






p. 


.5 


1 


0000 


b_total 


2.318431 


p. 


.6 


1 


0000 


Sigma_total = 


16.886438 a.u 


p. 


.7 


1 


0000 


Sigma_total = 


4.728683e-16 



According to the equations (fTTT i and (fTZt the A^*^ parameter is 

used only if v < Uy. In this case the A^^^ values are printed out 
in the Nef f column. Otherwise, when v > Uy, only a dashed 
line is printed. 

For calculations of the deposited energy T{b) the program 
prints out the b-mesh information and two columns (/j-values 
and T(b)-values). 

b_grid: [0.000; +0.010; 3.000] Npoints = 301 

Deposited energy T(b) calculation 
b T(b) 
0.000000 706.586257 
0.010000 706.333016 
0.020000 705.566015 



For the m-fold cross-section cr„, calculations, it prints out 
input values of ionization potentials, /j-grid data for the inte- 
gral (l22l i and values of the deposited energy. Once T{b) energy 
is calculated, the cr,„ data are evaluated and printed out. To 
check out the calculation the sum rule (|20] | is applied to the cr„, 
values and printed after the partial cross sections. 

m-fold cross-sections calculation 
I_l = 1.2660 a.u. = 34.450 eV 

I_2 = 1.7787 a.u. = 48.400 eV 



2.980000 0.216640 
2.990000 0.210448 
3.000000 0.204421 
Output T(b) into the file: 
energy _Tb_ba2+_o_out .txt 



I_28 = 62.0329 a.u. = 
I_29 = 65.7813 a.u. = 
I_30 = 69.5297 a.u. = 

m-fold Probability grid: 
N_points = 232 



1688.000 eV 
1790.000 eV 
1892.000 eV 



For the cr,£„ calculations the program prints out the input 
parameters, bisection search history (|39^ . and interpolation 
data ( l40l i. Finally, it prints out the /j^ax and cr,^, values. 

Total cross-section calculation 

II = 1.2660 a.u. 

b_min = 0.0000 b_max = 11.0000 

Bisection search: 



b_max 

grid_step = 
Save points: 
b 

0.0000 
0.0100 
0.0200 
0.0300 



2.3184 
0.0100 

T(b) 
706.5863 
706.3334 
705.5674 
704.2550 



2.2885 
2.2984 
2.3084 
2.3184 



1.3606 
1 . 3284 
1.2969 
1.2660 





b 


T(b) 


- 11 


m-fold 


Cross-sections : 





000000 


705. 


320243 


m 




a.u. 




cm2 


11 


000000 


-1. 


266014 


1 


6 


512001 


1 


823545e-16 


5 


500000 


-1. 


265965 


2 


3 


688431 


1 


032866e-16 


2 


750000 


-0. 


851313 


3 


1 


949328 


5 


458672e-17 


1 


375000 


8. 


266428 


4 


1 


146097 


3 


209396e-17 


2 


062500 


1. 


021403 


5 





750239 


2 


100883e-17 


2 


406250 


-0. 


245276 












2 


234375 


0. 


280613 


26 





000030 


8 


382028e-22 


2 


320312 


-0. 


005746 


27 





000000 


2 


943365e-25 


2 


277344 


0. 


131173 


28 





000000 


1 


303180e-30 


2 


298828 


0. 


061201 


29 





000000 


2 


553841e-40 


2 


309570 


0. 


027356 


30 





000000 





OOOOOOe+00 


2 


314941 


0. 


010713 


sum 


16 


886438 


4 


728683e-16 
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Output Pm(b) into the file: 
probability_Pm_ba2+_o_out .txt 

Start time: Sun May 20 14:18:36 2012 
Stop time: Sun May 20 14:27:55 2012 

Example of the calculated m-fold cross sections versus colli- 
sion energies is given in Figure |3] 

6. Conclusions 

A description of the DEPOSIT computer code, based on the 
energy-deposition model, is presented. This approach gives an 
agreement with experimental data within a factor of two. Nu- 
merical implementation of the model is given as well. The al- 
gorithms realized in the program allows one to simulate ion- 
atomic collisions at low and intermediate energies for arbitrary 
ion-atomic colliding systems. The code is parallelized and can 
be run in single and parallel modes (on a user PC or supercom- 
puter cluster). The computational speed is proportional to the 
number of cores. 
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Ba + O at 2.5 MeV/u (v = 10 a.u.) 




0.5 1 1.5 2 2.5 

b[ao] 

Figure 1: Energy deposition Tip) (in a.u.) to Ba^^ ion by an O atom at a collision energy E = 2.5 MeV/u (v = 10 a.u.) as a function of the impact parameter b, 
calculated with the DEPOSIT code. Energies deposited onto different electronic (Slater) shells of Ba^* are shown together with the total energy deposition. The 
horizontal dashed hne h = 1.266 a.u. (34.45 eV) corresponds to the first ionization potential of Ba^^ and shows the minimal energy deposition required for the 
ionization process. The notation for the Slater shells, for example, 2(spf means the electron configuration of the Is^lp^ shell. 
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Ba + O at 2.5 MeV/u (v = 10 a.u.) 




0.5 1 1.5 2 

b[ao] 



Figure 2: Calculated probabilities P,n{b) for the m-fold ionization of Ba ions colliding with O atoms at energy E = 2.5 MeV/u (v = 10 a.u.) as a function of the 
impact parameter b, eqs. )23t - )25t . 
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0.1 
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Figure 3: The total and m-fold electron-loss cross sections of Ba-* ions colliding with O atoms as a function of ion energy for ejection of m = 1, 2, 3, 4 electrons 
calculated by the DEPOSIT code (solid symbols). Experimental data from 1531 are available only for the m-fold cross sections (open symbols). At energies 
E < 2 keV/u only one-electron losses occur, that is why total and one-electron loss cross sections are the same in this region. 
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